Metastable Resistance Anisotropy Orientation of Two-Dimensional Electrons in High 

Landau Levels 



o 
o 

(N 

> 
O 



K. B. Cooper 1 , J. P. Eisenstein 1 , L. N. Pfeiffer 2 , and K. W. West 2 
1 California Institute of Technology, Pasadena CA 91125 
2 Bell Laboratories, Lucent Technologies, 
Murray Hill, NJ 07974 



In half-filled high Landau levels, two-dimensional electron systems possess collective phases which 
exhibit a strongly anisotropic resistivity tensor. A weak, but as yet unknown, rotational symmetry- 
breaking potential native to the host semiconductor structure is necessary to orient these phases in 
macroscopic samples. Making use of the known external symmetry-breaking effect of an in-plane 
magnetic field, we find that the native potential can have two orthogonal local minima. It is possible 
to initialize the system in the higher minimum and then observe its relaxation toward equilibrium. 
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There is now strong evidence for the existence of a 
broad new class of collective phases of two-dimensional 
electron systems (2DESs) which resemble charge density 
waves (CDWs) These new phases, which have only 
been observed in very high quality samples and at very 
low temperatures, occur when several of the discrete Lan- 
dau levels produced by a perpendicular magnetic field B± 
are occupied by electrons. Their transport characteristics 
allow them to be readily distinguished from the fractional 
quantized Hall states at very high magnetic fields and the 
weakly disordered Fermi liquid at zero field. Most strik- 
ing are those phases which develop around half-filling 
of the uppermost, or valence, Landau level (LL). Corre- 
sponding to LL filling factors v — 9/2, 11/2, 13/2, and so 
on, these phases are characterized by a huge anisotropy 
in the longitudinal resistance of the 2DES which develops 
rapidly on cooling below about 150 mK This im- 

pressive signature has been widely interpreted as strong 
evidence for the unidirectional, or striped, CDW pre- 
dicted via Hartree-Fock theory 0, 0- In this theory, 
and its many extensions 6], at half- filling of high LLs 
the exchange-induced softening of the Coulomb interac- 
tion leads to phase separation of the 2DES into stripes 
in which the valence LL is alternately filled and empty. 
Transport anisotropy is expected since it is almost cer- 
tainly easier for charge to move parallel to the stripes 
than perpendicular to them. The period of the stripe 
pattern is estimated to be a few times the classical cy- 
clotron radius R c . 

A key aspect of the anisotropic phases which lacks un- 
derstanding is their consistent orientation with respect 
to the crystalline axes of the host GaAs material. In vir- 
tually all samples this orientation is such that the high 
or "hard" resistance direction is parallel to (110) while 
the low or "easy" direction is parallel to (110). Numerous 
suggestions of weak native rotational symmetry-breaking 
potentials have been advanced, but so far none has been 
clearly shown to be relevant. Early experiments demon- 
strated that a magnetic field component B» in the plane 
of the 2DES provides an important tool for studying 



this potential. Properly directed, an in-plane field can 
overcome the native symmetry-breaking potential and rc- 




A good understand- 
and, via comparison 



orient the resistive anisotropy 
ing of this effect now exists 
with the experiments, has shown that the strength of 
the native symmetry-breaking potential is only about 1 
mK/electron Here we present strong evidence that 
the native potential can in fact possess two distinct local 
minima, one favoring the hard resistance direction along 
(110) and one favoring it along (110). Which is the global 
minimum depends sensitively upon filling factor and in- 
plane magnetic field. We find it possible to prepare the 
system in a metastable state of the "wrong" orientation 
and observe its slow relaxation toward the equilibrium 
state determined via a field-cooling technique. 

Figure la shows the longitudinal resistance vs. perpen- 
dicular magnetic field B±_ of a high mobility 2DES at T 
= 100 mK. The 2DES resides in a 30 nm GaAs quantum 
well and has a density and low temperature mobility of 
N s = 3x 10 11 cm -2 and /i = 3 x 10 7 cm 2 /Vs, respec- 
tively. Resistance measurements are made, via conven- 
tional lock-in techniques, using ohmic contacts positioned 
at the corners and side midpoints of the 5 mm square 
sample. For the solid trace the average current flow is 
along (110) and the resulting resistance is denoted R xx , 
while for the dotted trace the current flow is along (110) 
and the resistance is R yy - (These conventions will be fol- 
lowed throughout this paper.) The data shown in Fig. 
la cover LL filling factors v = hN s /eB± between v = 4 
and v = 8, corresponding to the Fermi level lying in the 
spin-split third (N=2) and fourth (N=3) orbital LLs. In 
spite of the relatively high temperature at which these 
data were obtained, there is much evidence of collective 
phenomena in Fig. la. For example, near B±_ w 2.6 
and 2.9T both R xx and R yy vanish. Although these fea- 
tures (indicated by upward arrows in the figure) occur at 
fractional filling factors (y « 4.72 and 4.28, respectively) 
they in fact reflect integer QHE states 0, E3|- Distinct 
from the conventional integer QHE states at v — 4 and 
5, these new states are widely believed to be pinned tri- 
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Perpendicular Magnetic Field (T) 

FIG. 1: a) Longitudinal resistances R xx (solid) and R yy (dot- 
ted) at T — 100 mK. b) and c) Sweep direction dependences 
of R xx and R yy at T = 50 mK around v — 9/2. d) and e) 
Sample tilted by 6 — 7°. R xx and R yy show a strong and 
non-hysteretic anisotropy at v = 9/2 at T = 50 and 100 mK. 



angular "bubble" CDWs closely related to the unidirec- 
tional stripe phases at half filling. 

More importantly, around filling factors v = 11/2 in 
the N=2 LL and v = 13/2 and 15/2 in the N=3 LL, the 
longitudinal resistance exhibits the strong anisotropy as- 
sociated with the putative stripe phases, with (110), as 
usual, being the hard direction. Interestingly, however, 
near v = 9/2, where anisotropy is usually strongest, the 
resistances R xx and R yy are almost identical. As we will 
show, this observation does not imply a lack of local 
stripe order, but instead indicates that here the native 
symmetry-breaking potential has two nearly degenerate 
minima along the (110) and (110) directions. 

Figures lb and lc display the resistances R xx and 
Ryy, respectively, in the region around v = 9/2 at the 
lower temperature of 50 mK. Each panel contains two 
traces; one taken with the magnetic field sweeping up 
at dB/dt = 0.01 T/min and one with the field sweeping 
down at the same rate. Unlike the situation at 100 mK, 
both R xx and R yy exhibit a strong dependence upon the 
direction of the field sweep at 50 mK. This dependence is 
confined to the region around half-filling; in the flanks of 
the LL, including the re-entrant integer QHE states, no 
sweep direction dependence is seen. These data make it 



clear that hard and easy transport directions cannot be 
assigned at v — 9/2: On sweeping the field up the hard 
direction appears to be along (110); on sweeping down it 
seems to be (110). A similar effect is seen in the region 
around v = 11/2, even though at 100 mK (110) is clearly 
the hard axis. 

Figures Id shows that the application of a small in- 
plane magnetic field completely removes the sweep di- 
rection dependence of the resistance around v = 9/2. 
To obtain these data the sample was tilted by 9 = 7° 
around the (110) direction; at v = 9/2 this produces an 
in-plane field of B\\ — 0.34 T directed along (110). In this 
situation the hard axis of the resistive anisotropy is un- 
ambiguously (HO), in agreement with earlier tilted field 
measurements Fig. le demonstrates that the small 

applied in-plane field is sufficient to bring out a strong 
anisotropy at v — 9/2 even at 100 mK. 

Similar observations have been reported by Zhu, et 
al. 0]. I n their experiments the 2DES density was 
swept at fixed magnetic field. A strong dependence of 
the longitudinal resistance on the direction of the den- 
sity sweep was observed at v = 9/2 if the density was 
near N s — 2.9 x 10 11 cm -2 . For densities well below this 
value the hard direction at v = 9/2 was found to be 
(110), while for higher densities it was (110). Zhu, et al. 
did not propose an explanation for this density-driven re- 
orientation of the resistive anisotropy and suggested that 
their results could not rule out either a unidirectional or 
a bidirectional native symmetry-breaking potential. 

Our experiments, which employ a sample with a fixed 
density close the critical one identified by Zhu, et al. ^3] , 
strongly point to the existence of a bidirectional poten- 
tial. The sweep direction dependence of the resistivity 
shown in Figs, lb and lc proves that over a range of 
magnetic fields about v = 9/2 there are two orthogo- 
nal orientations simultaneously possible for the resistive 
anisotropy These two orientations likely correspond 
to two distinct minima of the symmetry-breaking poten- 
tial. Apparently the relative depth of these two minima 
is a sensitive function of filling factor. For magnetic fields 
slightly below v = 9/2, the data in Figs, lb and lc imply 
the potential minimum for stripes with hard axis along 
(110) is lower in energy. Conversely, for fields slightly 
above v = 9/2, the potential favors the hard axis along 
(110). If the time required for stripe domains to relax 
from one minimum to the other is sufficiently long, then 
this picture would naturally explain the sweep direction 
dependence of the resistances. The elimination of this 
hysteresis by an in-plane magnetic field is a consequence 
of the in-plane field's own symmetry-breaking effect fa- 
voring one of the two minima of the native potential. 

In order to make this model more convincing we now 
turn to the behavior of the resistance near v = 13/2. 
Here we have directly observed the slow relaxation of 
metastable anisotropy orientations toward an equilibrium 
configuration. Moreover, a field-cooling technique allows 
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FIG. 2: a) and b) Resistances R xx (solid) and R yy (dotted) 
near v — 13/2 at T = 50 mK showing extreme sweep direction 
dependence. Sample tilted by 6 — 7°. c) Time evolution of 
R xx and R yy following field sweep down to B±— 1.89 T at 
50 mK. d) R xx (solid dots) and R yy (open dots) after cooling 
from 150 to 50 mK at fixed magnetic field, e) Initial time 
dependences of R xx following field sweep down to B±=1.89 
T at 50 mK. Dashed lines are field-cooled values of R xx . 



for unambiguous identification of the equibrium configu- 
ration, avoiding the complications of hysteresis and slow 
relaxation times. 

The data in Fig. la show that at T = 100 mK and 
Bn — a large resistive anisotropy is present at v = 13/2, 
with the hard direction along (110). On cooling to 50 mK 
we find that a slight magnetic field sweep direction de- 
pendence does develop, but the hard transport direction 
remains (110) regardless of the sweep direction. Within 
our model of a bidirectional symmetry-breaking poten- 
tial, this implies that while the relative depth of the two 
local minima varies with magnetic field around v = 13/2, 
which minimum is lowest in energy docs not change. The 
small hysteresis observed presumably arises from a small 
non-equilibrium population of stripe domains with the 
"wrong" orientation. If this is the case, the application 
of an appropriate in-plane magnetic field ought to bring 
the two minima closer together and exacerbate the sweep 
direction dependence. As Figs. 2a and 2b show, this is 
precisely what we observe. With the sample once again 
tilted to 6 = 7°, producing B\\ = 0.23 T along (110) at 
v = 13/2, the observed sweep dependence of the resistive 



anisotropy is extreme. A large anisotropy of very simi- 
lar aspect is seen in both Fig. 2a and 2b, only the hard 
and easy directions are interchanged. As with v = 9/2, 
sweeping the magnetic field up favors the hard direction 
along (110), while sweeping it down favors (110). 

If the sweep direction dependence of the resistive 
anisotropy is due to striped domains getting trapped in 
the "wrong" direction as the magnetic field is swept into 
the region of half-filling, owing to a filling factor depen- 
dence of the relative depth of the two potential minima, 
then it is reasonable to expect to observe temporal relax- 
ation of the resistance, provided that the energy barrier 
between the two wells is not too large. Figure 2c shows 
that just such relaxation is observed. In the figure R xx 
and Ryy at T = 50 mK are plotted vs. time, over a 19 
hour interval. These data were obtained by first sweep- 
ing the magnetic field down to B± = 1.89 T, stopping, 
and then recording the subsequent time evolution. This 
stopping field was chosen because it is in the center of 
the anisotropic region about v = 13/2. Immediately after 
stopping at 1.89 T, R yy « 800 fi, while R xx is only about 
s=s 7 Vt. Subsequently R yy falls quickly at first and then 
slows down. Meanwhile, R xx slowly rises. After about 
11 hours the two resistances cross and (110) becomes the 
hard direction. As time evolves further, the gap between 
R xx and R yy grows steadily. These data clearly suggest 
that the initial anisotropy orientation present when the 
field sweep was stopped (hard axis along (110)) was a 
metastable one and that the equilibrium orientation has 
the hard axis along (110). 

To confirm this scenario we have employed a field- 
cooling technique to determine the equilibrium config- 
uration of the anisotropy. In this technique the magnetic 
field is swept into the v = 13/2 vicinity at T = 150 mK. 
At this high temperature no anisotropy is evident in the 
resistances and no hysteresis is observed. Once the field 
sweep is stopped, the sample is allowed to cool to T = 50 
mK and both R xx and R yy are then recorded. Figure 2d 
shows the results of this procedure. The solid dots rep- 
resent the T = 50 mK values of R xx ; the open dots R yy - 
We emphasize that no time dependence of these values 
is observed after cooling to 50 mK and no dependence 
upon the direction of the 150 mK field sweep is found. 
It is clear from Fig. 2d that for most of the magnetic 
field range around v — 13/2, the equilibrium hard direc- 
tion is (110). Only in a very narrow range on the high 
field side of the anisotropic region does (110) appear to 
assume that role. 

That the field-cooled resistances shown in Fig. 2d are 
the equilibrium ones is demonstrated in Fig. 2e. Each 
panel displays both the time dependence of R xx over a 10 
minute interval following a downward field sweep to B±_ = 
1.89 T at T = 50 mK, and the field-cooled value (dashed 
horizontal lines) of the same resistance obtained via the 
field-cooling technique described above. The data show 
that the R xx values observed after the 50 mK field sweeps 



4 



Field-cooled to 50mK 




2.6 2.8 3.0 



B ± (T) 

FIG. 3: Rxx (solid dots) and R yy (open dots) in vicinity of 
v = 9/2 after field-cooling to T — 50 mK. Sample is untilted. 

do indeed relax toward the equilibrium values found after 
cooling from 150 mK at fixed field. Not surprisingly, the 
relaxation time is strongly temperature dependent: ft is 
many hours at 50 mK, several minutes at 70 mK, and 
almost immediate at 90 mK. Assuming the relaxation 
rate is determined by thermal activation over the barrier 
between the two potential minima, a rough estimate of 
the barrier height is C4 ~ 1 K. 

Applying the field-cooling technique to the resistances 
around v = 9/2 in the N = 2 LL yields the data in Fig. 
3. These 50 mK data were obtained with J5| | = and 
therefore only the native symmetry-breaking potential is 
present. These equilibrium results further support our 
model of a bidirectional native symmetry-breaking po- 
tential. Only a modest level of resistive anisotropy is 
present in the window 4.35 < v < 4.65 centered on half- 
filling. For 4.35 < v < 4.57, R yy exceeds R xx , while 
for 4.57 < v < 4.65 R xx exceeds R yy . These results 
suggest that the two orthogonal minima of the native 
potential are nearly degenerate throughout the 9/2 re- 
gion, a conclusion consistent with the lack of anisotropy 
at T = 100 mK evident in Fig. la. Moreover, the 
hard resistance direction being along (110) for magnetic 
fields above v = 4.57 and along (110) for fields below 
this point, makes sense of the hysteresis data shown in 
Figs, lb and lc. Sweeping into the 9/2 region from 
high fields would initially generate striped domains with 
hard axes along (110), while approaching from low fields 
would initially produce domains with hard axes along 
(110). Given the demonstrated long relaxation times of 
metastable anisotropy orientations at low temperatures, 
hysteresis is obviously expected. 

The existence of two orthogonal local minima in the na- 
tive anisotropy potential, one favoring the hard axis along 
(110) and one favoring (110) implies that in addition to a 
term in the native potential proportional to cos(20) (with 
being the angle between the hard transport direction 
and (110)), there is also a term proportional to cos(40). 



The sign of the cos(20) term determines whether (110) 
or (110) is the preferred hard axis direction. Prior ex- 
periments 0| suggest that the sign of this term changes 
at a critical density. Our results are consistent with this 
but reveal additional dependences on LL index and the 
precise filling within each spin-resolved LL. The physical 
origin of the cos (20) term remains unknown (llj . Unlike 
the cos (20) term, a cos (40) term is not prohibited by the 
crystal symmetry of bulk GaAs. Our results demonstrate 
that this term leads to dramatic metastability phenom- 
ena when the cos(20) term is small. The latter occurs 
when the density is near the critical one 0] or when an 
appropriate in-plane magnetic field is applied. One inter- 
esting recent suggestion 1S that the piezoelectricity 
of GaAs produces a cos (40) term with minima favoring 
(110) and (110). The barrier separating these degener- 
ate minima is roughly 1 mK/electron. Comparing this 
with our crude estimate J7& ~ 1 K for the height of the 
barrier which inhibits relaxation of mis-oriented stripe 
domains leads us to speculate that each domain contains 
about 10 3 electrons, corresponding to lateral dimensions 
of order Xfiia. 

In conclusion, we have presented evidence that the na- 
tive symmetry-breaking potential responsible for orient- 
ing the anisotropic phases near half-filling of high Landau 
levels has two orthogonal local minima, one favoring the 
hard resistance direction along (110) and one favoring 
(110). We have shown that it is possible to witness the 
dynamic relaxation of the anisotropy from metastable 
states to an equilibrium orientation separately deter- 
mined via a field-cooling technique. 
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